PACS numbers: 47.50.-d, 47.20.Ft, 47.54.-r 

Thermodynamic potential of a mechanical constitutive model for 

two-phase band flow 

Katsuhiko Sato^, Xue-Feng Yuan^ and Toshihiro Kawakatsu^ 
^ Department of Physics, Tohoku University, Sendai 980-8578, Japan and 
^ School of Chemical Engineering and Analytical Science, 
The University of Manchester, Manchester, M60 IQD, United Kingdom 

(Dated: November 10, 2008) 

Abstract 

Starting from a simple mechanical constitutive model (the non-local diffusive Johnson-Segalman 
model; DJS model), we provide a rigorous theoretical explanation as to why a unique value of the 
stress plateau of a highly sheared viscoelastic fluid is stably realized. The present analysis is based 
on a reduction theory of the degrees of freedom of the model equation in the neighborhood of a 
critical point, which leads to a time-evolution equation that is equivalent to those for first-order 
phase transitions. 
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Highly sheared viscoelastic fluids are known to have spatially inhomogeneous structures, 
i.e., "shear bands" [l, 2], where the fluid is separated into two regions with different values 
of the velocity gradient. In this phenomenon, a stress plateau is observed in the stress-strain 
curve (SS curve). 

This phenomenon has been explained using an N-shaped SS curve, as depicted in Fig. 
1 [sl. In the region of the shear rate where the SS curve has a negative slope (referred 
to hereafter as the negative-slope region), the uniform flow is unstable, and then the fluid 
separates into two stable domains having two unique values of shear rates, denoted by 71 
and 72 (71 < 72). As the applied shear rate 7 increases, the relative volume fraction of the 
domain with 72 increases. This scenario enables us to explain the appearance of the shear 
stress plateau using a consideration similar to thermodynamics. The validity of this scenario 



y,a. 



has been generally conflrmed by direct experimental observations 

On the other hand, there are some theoretical disadvantages associated with this sce- 
nario. One such disadvantage is the fact that the origin of the unique value of the stress 
plateau, Us, is unknown. (In experiments, the value of the stress plateau is known to be 
uniquely determined independently of the flow history.) A solution to this problem is given 
w adding a non-local term (a diffusion term) to the constitutive model. In References 
3, y, the uniqueness of the value of (Xs is demonstrated numerically, and in Reference 
10l | a selection rule for a^, which is analogous to the Maxwell equal area construction, is 
provided analytically. 

Although the problem of the unique determination of as was proven in the above refer- 
ences, another important problem remains, i.e., the global stability of the banded state. As 
stated above, in the negative-slope region the banded state is flnally realized because of the 
instability of the homogeneous flow. However, in the remaining regions, (a)-(c) and (d)-(b) 
as indicated in Fig. 1, the SS curve has a positive slope and so the homogeneous flow is 
stable (at least locally). Thus, it is not guaranteed that the banded state is flnally realized 
in these regions. In other words, we cannot determine whether the homogeneous flow or the 
banded flow is more stable. This is because we have not yet had an evaluation function such 
as a thermodynamic potential for the models describing the shear banding. 

In the present article, we will demonstrate the metastability of the above-mentioned ho- 
mogeneous flow by deriving an evaluation function from a simple mechanical model, called 
Johnson-Segalman (JS) equation with a diffusion term (DJS equation). The DJS model 



2 



Danding, with which characteristic ex- 
Here, the term "mechanical" 



is a widely accepted model paradigm of shear 
perimental results are well reproduced numerically 
indicates that the model is constructed in a purely mechanical manner, i.e., the model does 
not have an explicit thermodynamic potential. While there are some other extended models 
on shear bands 



Q, Q, as a first trial on this issue we will adopt the simplest DJS 



[HI, 

equation. In fact, the method we use in the present study, which is based on the reduction 
theory developed in the field of non-linear dynamics |20i], is also applicable to a wide class 
of models. 

We start with the following governing equations for a viscoelastic fiuid: 

V-v=0, (1) 
p{^^+v-V)v=V-{^ + 27]D~pI), (2) 

(— + V ■ V)'E=a(D ■ S + S ■ i:)) + (n ■ S - S • n) + 2GD 2DoV^D. (3) 

at T 

Equation ([1]) represents the incompressibility condition of the fiuid, where v is the fiuid 
velocity. This condition holds in usual fiuids and guarantees that the density p is constant 
throughout the fiuid. Equation ([2]) represents the momentum balance of the fiuid, where 
the total stress tensor denoted by cr is assumed to be composed of two contributions as 
cr = S + 21] D, where SI is the contribution from the polymeric components and 2riD is 
the contribution from the solvent. We assume that the polymeric component S obeys the 
DJS equation given by Eq. Q, while the solvent is assumed to be a Newtonian fiuid so 
that its contribution is a product of the viscosity rj and the symmetric part of the velocity 
gradient tensor D = ((Vv)^ + Vv)/2. The quantity p is the pressure, and / is the unit 
tensor. Equation is the DJS equation, where a is a parameter that satisfies < a < 1. 
This parameter a is called the "slip parameter", which represents the degree of non-affine 
deformation under the fiow [l4]. O is the antisymmetric part of the velocity gradient tensor, 
defined by O = ((Vf)"^ — Vv)/2. The coefficients G, t and Dq in Eq. ([3]) are scalar 
constants having the dimensions of stress, time, and (stress x (length)^), respectively. The 
last term V^D is the diffusive stress term, the components of which are '^^d'^Djk/dxf for 



any j and k. In References js, 9|, the non-local term of the DJS equation is assumed to have 
the form V^S. The physical meaning of this form can be understood as a diffusion of the 
stress tensor originating from the diffusive flux of the polymeric components. Although the 
physical meaning of V^-D appears to be less obvious, it will be shown that the flnal results 



of the present analysis are unchanged both of these choices. (See the last paragraph of this 
article.) 

Using this setup, we consider a one- dimensional simple shear flow, i.e., a flow contained 
between two parallel plates moving with a constant relative speed Vq along the x-axis. The 
velocity distribution of the fluid v = {vx,Vy,Vz) in such a situation is given as = Vx{y,t) 
Vy = Vz = with the boundary conditions Vx{0,t) = and Vx{L,t) = vq for any time t, 
where L is the distance between the two plates. For the solvability of the model, we need 
additional boundary conditions. We require d'^Vx{y,t) /dy'^\y=o^L = for any time t. This 
boundary condition corresponds to the situation in which the polymeric components are 
impenetrable at the walls of the plates Q. 

To simplify the equations, we use the same non-dimensionalization procedure introduced 
in early investigations j^, [ij], i.e., t/r t, y/L y, pLi^/r'^G p, rj/Gr —>■ rj and 
Dq/2GL^ Di. By applying this non-dimensionalization, we obtain three non-dimensional 
quantities, K{y,t), N{y,t) and S{y,t), defined hy k = \J\ — o?t^^, N = [(1 — a)T.xx — 
(1 + a)Tjyy]/2G, and S = Vl — o?^xy/G- Then, Eqs. ([I])-© reduce to the following set of 
equations for N and S: 

f-.S-A'. (5) 

. Here, Eq. (|4]) is obtained by differentiating the x-component of both sides of Eq. ([2]) and 
eliminating Vx- In this non-dimensional system, the boundary conditions are given as 



-.-kN + k-S-D, — k. (6) 



/ K{y, t)dy = K and 
Jo 



= 0, (7) 

j/=0,l 



dy 

where k is a constant given hj R = t\J\ — o^vq/ L. In the following analysis, we assume 
that the quantity Di is an extremely small positive number compared to unity, which may 
be justified by the fact that Di oc l/L^. This assumption corresponds to the situation in 
which the width of the interface between the shear bands is much smaller than the system 
size because in this model the width of the interface is given by a quantity proportional to 

A trivial y-independent stationary solution of Eqs. (jl])-© that satisfies the boundary 
conditions of Eqs. ^ \s n = N = R'^ /{I + k^), and S = R/{1 + R"^). This solution gives 



a trivial value of the total stationary shear stress a = ^-^_^-2-^ + riR. When 77 > 1/8, a is a 
monotonic function of K, and when 77 < 1/8, a is a non-monotonic N-shaped function of k 
that has a region with a negative slope. The negative-slope region of the DJS equation is 
given by k„„ < R < Hn+ with nn± = a^/— 1 + 4= ± "^^"3^- '^^^ trivial stationary solutions 
have been proven to be unstable in this region When rj approaches 1/8, the negative- 

slope region converges to a single point R = ^/S, which allows us to define the 

"critical point" of the DJS equation as {R*,r]*) = (y^, 1/8). At this critical point, the 
conditions da/dR\fi=fi* = and (Pa / dR'^\f^=f^* = hold simultaneously. 

A non-trivial stationary solution of Eqs.(j4])-(l6l) in the negative-slope region for rj < 1/8 
is obtained in the case with the V^-D term in Eq. ([3]). By assuming a stationary solution 
of Eqs.(!4])-(l6j) under the boundary condition f (y, t)|j;=o,L = const, we obtain 

k{1 - K{as - r]K)) - {as - rjK) - Di-^ = 0. (8) 

We can identify this equation with an equation of motion for a particle with position k at 
time y moving in a quartic potential. This equation has a non-trivial solution that satisfies 
the boundary conditions Eq.Q only when 



crs = 3^-v' + l (9) 

Note that this value of does not depend on R. At this value of ag, the non-trivial solution 
of Eq. ([H]) is given as 

K4- — K n 



i^{y) = i^oiy) = ^ ^ tanh 



y 



(10) 



, where k± = — — ^ '' — . This non-uniform solution describes a shear banding, in which 
the relative volume fraction of the regions with lower shear rate is given by (k^ — R) / — k_) 



and the width of the interface between different regions is given hy ^ = ^y Di/(1 — 8r]). From 
the symmetry of Eq.(l8]), we know that K,{y) = Ko(1 — y) is also a non-uniform solution of Eq. 
(IH]). It is worth noting that the analytical results of Eqs.Q and ffTOl) are in good agreement 
with the numerical results for the same model [?[. 

As mentioned in the introduction, a linear stability analysis shows that the trivial homo- 
geneous solution is unstable in the negative-slope region < k < k„+), and therefore a 
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nontrivial solution should emerge. However, in the remaining regions < k, < Kn- and 
K,n+ < n < where a has a positive slope, there is no guarantee that the non-trivial solu- 
tion is finally realized because the homogeneous fiow is still stable (at least locally) in these 
regions. Therefore, for a complete understanding of the stress plateau, it is not sufficient to 
determine the value of ag uniquely and obtain the corresponding non-trivial solution. The 
metastability of the homogeneous fiow in the positive-slope regions must be demonstrated. 
To do this, we shall reduce the DJS equation in the neighborhood of its critical point. 

The reduction procedure is as follows. First, we ignore the inertia term (the adiabatic 
approximation) in Eq. (j4]) because the Reynolds number is small for usual viscoelastic 
fiuid fiows. We then have an approximate expression for k, as K,{y, t) ~ k — e{S — Sq), where 
e = l/r], and 5*0 is the mean value of S over the spatial coordinate y, i.e., So{t) = S{y, t)dy. 
Substituting this approximate expression for k into Eqs.([5]) and ([6]) yields 

--{R-e{S-So))S-N, (11) 



dt 

_ e(^ - ^o))(l -N)-S + d^, (12) 



where d = Die. Next, in order to clarify which mode is dominant in the critical region, 
we carry out a linear stability analysis at the critical point {R*,e*) = (-\/3, 8). Linearizing 
Eqs.f fTT]) and f|T2l) with respect to the stationary solution at the critical point, {N*,S*) = 
(3/4, -\/3/4), in terms of the deviations sq, uq, s, and n defined by sq = Jq{S — S*)dy, 
riQ = Jq{N — N*)dy, s = S — S* — sq, and n = N — N* — uq, we have 

— - = -no + VSso -TTT = -f^- 

dso ds (1^) 

= -VSno - So { di ^ -V3n - 3s. 

Here, we have neglected the diffusion term because, according to the reduction theory for 
partial differential equations, the diffusion term should be regarded as a small perturbation 
to the uniform system j20|. The eigenvalues of the former set of equations are — 
and those for the latter set of equations are and —4. The corresponding eigenmodes of 
the latter set of equations are 



h-V3n + s), (14) 



2 



^=l^n + V3s), (15) 
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respectively. Thus, we find that, in the vicinity of the critical point, the variable 6 is the 



unique slow one of the system and the other variables no, sq and ip are solved by (p |2Q|, l21 |. 
Based on this observation, we first rewrite the original Eqs. lITT]) and (IT^ in terms of no, So, 
(p and ip as 



S*iR-R*)-no + Rso-^{i(p+V3ijy), (16) 



^^0 

^=(^ _ _ iv*) - Rn^ - so + |((-v^0 + ^)(0 + v^^)), (17) 

^=^(-v^/i + h) + ^5J(0 + v^V'), (18) 

%=\{h + + 4^9J(0 + 73^), (19) 

where the bracket denotes an integral over the entire domain of y, i.e., (■) = -dy. The 
quantities /i and /2 are functionals of no, So, 0, and ip, defined by /i = {— n + {—eS* + 
k)s — esos — e(s^ — (s^))}|p and /2 = {— Kn — (e(l — N*) + l)s + enQS + e{hs— (ns))}|p, where 
the subscript p indicates the substitutions n = |(— -\/30 + -i/^) and s = i(0 + VSip) in each 
expression. Note that these Eqs.(fT6l)-(fT9l) are still equivalent to the original Eqs.l fTTI) and 
(fT2l) . Next, we assume that the variables no, sq, and ip are functionals of (p and functions of y, 
denoted by no[0, y], so[(p, y] and ip[(p, y], and define a functional G[(p, y] as the right-hand side 
of Eq. ffTSl) . Substituting these forms into Eqs. f|T6l) - f|T9l) gives four algebraic equations for no. 
So, ip, and G with independent variables (p and y. Here, our task is reduced to determining 
the forms of these four functionals. Introducing two small parameters k and e defined as 
k = R — R* and e = e — e*, (|e| ^ 1, \k\ ^ 1) and realizing that the difference between 
the shear rates of each domain is on the order 0{^) (e > 0), which is known from the 
expanded form of a near the critical point, we can proceed to the perturbative calculation 
systematically. The results are no = ^k - ^ ~ 2(0^), ^o = -|, V = ^ - (0^ - (0^)) and 
G= (| — |A;^)0 + 3A;(0^— (0^))— 4(0^ — (0^)). Here we have used the property (0) = 0, which 
results from its definition (IT^ . In the above calculation, we have neglected the diffusion 
terms because of the smallness of d. If we retain the diffusion terms up to the lowest order 
in d, we have the following time evolution equation for 0: 

f = - + 3M0^ - (0^)) - 4(0^ - (0^)) + (20) 

This equation expresses the main result of the present research. The dynamic variable is 
essentially an "order parameter" that which measures the inhomogeneity of the viscoelastic 
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stress field and includes the contributions from both the shear stress and the first normal 
stress difference. Equation fl20l) has the same form as that of the time-dependent Ginzburg- 
Landau (TDGL) equation for first-order phase transitions in thermodynamics 22, |23|] . A 
non-trivial stationary solution of Eq. fl2U]) for e > is readily obtained, which is consistent 
with the exact solution of Eq. (11 01) to the lowest order in e. An important advantage of the 
reduced equation, Eq. (!20|) is that we are able to prove the metastability of the homogeneous 
flow, which is readily performed by conventional means, e.g., a common tangent construction. 
A detailed discussion on the properties of Eq. (!20!) will be reported in a future paper |24l |. 

Finally, we would like to comment on our choice of the diffusion term. In Eq. ([3]), even if 
we replace the diffusion term DqV^D with the usual form D2V^'S, the reduced TDGL-type 
equations is unchanged. The only difference is the diffusion constant appearing in fl2Ul) . 
where d/4 is replaced by tD2/L'^. 



Acknowledgments 

The authors would like to thank T. Shibata, Y. Kuramoto, T.Yamaguchi, T. Chawanya, 
N. Uchida, T. Ohta, H. Watanabe, and T. Koga for their helpful discussions and com- 
ments. The present study was supported by Grant-in- Aid for Scientific Research on Priority 
Area "Soft Matter Physics" from the Ministry of Education, Culture, Sports, Science, and 
Technology of Japan. 



8 




FIG. 1: Figure 1: Schematic representation of shear banding. In the regions of shear rate where 
the SS curve has a positive slope, 71 < 7 < 73, and 74 < 7 < 72, both the homogeneous flow and 
banded flow are stable. Hence, in order to determine which state is finally realized in these regions, 
it is necessary to have some evaluation function to compare the stabilities of these states. 
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